Autangle: A case of Quantum Narcissism? 
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In this paper we ask a common psychological question and provide a physics answer: "Looking 
into a mirror can one get entangled with one's image?" This is not a frivolous question; rather, it 
bears on the effect of boundaries on the behavior of quantum entanglement between a harmonic 
oscillator and a quantum field, a basic problem of interest in proposed mirror-field superposition and 
related experiments in macroscopic quantum phenomena, as well as atomic fluctuation forces near a 
conducting surface. The object's internal degree of freedom is modeled by a harmonic oscillator and 
the presence of a perfectly reflecting mirror enforces the Dirichlet boundary conditions on the quan- 
tum field, restricting the latter to a half space. By assuming a bilinear oscillator-field interaction, 
we derive a coupled set of equations for the oscillator's and the field's Heisenberg operators. The 
former can be cast in the form of a quantum Langevin equation, where the dissipation and noise 
kernels respectively correspond to the retarded and Hadamard functions of the free quantum field 
in half space. We use the linear entropy as measures of entanglement between the oscillator and 
the quantum field under mirror reflection, then solve the early-time oscillator-field entanglement 
dynamics and compare it with that between two inertial oscillators in free space. At late times 
when the combined system is in a stationary state, we obtain exact expressions for the oscillator's 
covariance matrix and show that the oscillator-field entanglement decreases as the oscillator moves 
closer to the mirror. We explain this behavior qualitatively with the help of a mirror image and 
provide an answer to the question raised above. We also compare this situation with the case of two 
real oscillators and explain the differences. 



* zhour@umd.edu 
t rbehuninOlanl.gov 
t sylm@cc.ncue.edu.tw 
§ blhu@umd.edu 



2 



I. INTRODUCTION 

In this note we ask a common psychological question and provide a physics answer: "Looking into a 
mirror can one get entangled with one's image"? According to storybooks the answer for the evil queen 
Q seems to be yes, a case of typical royal narcissism. We want to find out the answer from physical 
considerations: here entanglement refers to quantum entanglement, and autangle means self-entanglement 
, here referring to the entanglement of a real physical object Q with its mirror image. We hasten to 
add that since the image is not a physical object one cannot define, let alone calculate, the entanglement 
between a physical object and an unphysical construct. However, there exists entanglement between an 
oscillator and a quantum field, and we can ask how this is altered if a mirror is present. A perfectly 
reflecting mirror imposes Dirichlet boundary conditions on the field along the mirror surface, restricting its 
existence to a half space. These are well-defined problems, belonging to the broader inquiry into the effects 
of boundaries and topology on the quantum field and on quantum entanglement of objects, oscillators and 
mirrors, coexisting with and/or mediated by the quantum field. In fact, they bear on issues of importance to 
quantum information and macroscopic quantum phenomena. Amongst research problems of current interest, 
we mention: 1) Entanglement between two two- level (2L A) atoms via a common quantum field already shows 
diverse quantum entanglement dynamics behavior |2|, ranging from sudden death, touch of death, revival 
to staying always alive, and features such as dynamical generation, protection, and transfer of entanglement 
between subsystems. Model studies of oscillator-field entanglement with experimental ventures have been 
carried out [3j as well as oscillator-field entanglement in conjunction with measurements in LIGO detectors 
[H, the latter also serving as preparatory studies for mirror-field superposition [j| in macroscopic quantum 
phenomena. 2) Entanglement not only changes in time but it also depends on spatial separation as shown in 
model studies for 2LAs 0-Q and harmonic oscillators [9| interacting with a common quantum field. These 
results are of both theoretical interest in understanding what quantum nonlocality means, and practical 
value, such as for the design of quantum gates and applications to quantum teleportation. Extending these 
findings to two and three dimensional systems would enable one to define quantum entanglement domains 
[l2l,[l3j and geometric effects of entanglement (l4lll5j. The former refers to the effective domain on a surface 
or body where entanglement is induced by the presence of an oscillator, not unlike the induced surface charge 
density on a conductor or dielectric plate. The latter refers to specific geometric patterns in the arrangements 
of oscillators where their entanglement strength can be maximized. 

Now, returning to the problem under study, one may wonder, is it really necessary to carry out such a 
serious calculation in order to answer a simple query posed at the beginning: Is there entanglement between 
the queen and her image? The answer is yes, both to the query and to the necessity of a bona fide calculation. 
As we cautioned ab initio, an image is not a physical object. Thus we try to avoid invoking an image in 
this calculation so that no unphysical assumptions or intuition are brought in unnecessarily or unknowingly 
in the derivations. However, at the end, by inspecting the results from the oscillator-field (under mirror 
reflection) entanglement we see that one can use the notion of an image to describe the process, in fact, come 
up with a simple qualitative explanation of how it depends on distance. 

The paper is organized as follows: In Section. II we set up the problem and derive a formal set of equations 
for the Heisenberg operators of the oscillator's internal degrees of freedom and the quantum field. We impose 
the boundary conditions introduced by the presence of the mirror and explain how the commutation relations 
are altered. Next we calculate the retarded and Hadamard Green functions of the altered field configurations 
and derive a quantum Langevin equation for Q(t) including the back-reaction of the altered field. We seek 
solutions to this equation at late times and calculate the covariance matrix of the oscillator's canonical 
variables at late times when the combined system is in a stationary state. Details are contained in Appendix 
[Al In Sec. Ill we introduce the linear entropy and the von-Neumann entropy as measures of quantum 
entanglement and calculate the early-time dynamics of entanglement between the oscillator and the field in 
the oscillator-mirror setup described above. We obtain plots of how the entanglement between the oscillator 
and field in the half space with the mirror reflection evolves with time, depending on the distance between 
the oscillator and the mirror. Finally in the discussions section we explain in what sense can one describe 
this situation in terms of entanglement of the oscillator with its mirror image and why it decreases as the 
oscillator moves closer to the mirror, a somewhat counterintuitive finding. We also compare this situation 



We coin this term from two sources of inspiration: new words like 3-tangle and old ones like autism. A more learned word 
'ipso-tangle' was suggested by Prof. Brill. 
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with the case of two inertial oscillators in free space (calculations placed in Appendix [Bj and explain their 
physical differences. 



II. OSCILLATOR INTERACTING WITH FIELD UNDER MIRROR BOUNDARY CONDITION 

A. Dynamics of Oscillator-Field System with Mirror 

As described in the Introduction, the gist of the matter for this problem is to quantify the change in 
some entanglement measure between the oscillator and the field with and without the presence of a mirror. 
Technically our calculation involves two parts: 1) find the back-reaction of the field configuration altered by 
the mirror on the oscillator. 2) find the correlation functions of the resultant oscillator's canonical variables 
including the field's influences. Thus our target is the covariance matrix for the oscillator which incorporates 
all back reactions from the quantum field. We work with the Heisenberg equations of motion which efficiently 
provides the oscillator's field-influenced dynamics. We can then find the oscillator's late time behavior and 
quantify it's entanglement with the modified field. 

Consider an oscillator located at position Xq at a vertical distance of L/2 from the mirror plane at x 3 = 0. 
The total action of the system describing this oscillator with internal degree of freedom Q interacting linearly 
with a massless scalar field confined to the half space defined by X3 > is given by: 

S[Q,Q;$,d^] = ^M Q J dt (Q 2 -Sl 2 Q 2 ) + ^ J dt jf A 9,<D9"$ + A Q j dt Q{t)${x Q ,t). (2.1) 

Due to the linearity of our system the equations of motion for Heisenberg operators (carrying hats) have the 
same form as the equations for the corresponding classical variables: 

M Q Q{t) + M Q n 2 Q(t) =A Q f d 3 x $(x Q) i), (2.2) 

Jx 3 >0 

□$(x,i) =A Q( 5 3 (x- Xg)Q(f), and ®(x\\,x 3 = 0,t) = 0. (2.3) 

Equation (|2.3|l specifies the Dirichlet boundary (DBC) conditions imposed on the field where the mirror 
surface is located, namely, in the 2; 3 ~ plane and X| = {x\, £2)- To obtain the back-reaction of the field on 

the oscillator we first solve (|2.3|) and then plug the solution into (|2.2j) . The solution for $(x, i) is given by 

<l(x, f) = $ (x, t) + X Q J dt'G% t (t, x; t', x Q )Q(t') (2.4) 

where $o( x , t) is the homogeneous solution to equation (|2.3[) which describes the dynamics of the source- free 
field (without Q) in the presence of the mirror and G* et (x, t; y, t 1 ) is the retarded Green's function for the 
field. Plugging the solution for the field operator into the equation of motion for the oscillator we obtain the 
following equation governing the dynamics of the oscillator with the effects of the field already incorporated. 
This is what we mean by the 'field-influenced' dynamics of the oscillator 

M Q Q(t) + M Q n 2 Q(t) -\%J dt'G? et (t, x Q ; t\ x Q )Q(t') = A Q <l>o(x Q , f). (2.5) 

As will be explained in further detail later the term containing the retarded Green's function describes 
how the oscillator transfers energy to the field and the right hand side acts similarly to a Langevin forcing 
term describing how quantum field fluctuations drive the oscillator. 



B. Field Modified by Mirror and Propagators in Half-space 

In order to solve the equation of motion for the oscillator (|2.5p we need to know the field operators defined 
in the half space. Following this procedure we first write down the free field operator $0 which satisfies the 
Klein-Gordon equation and vanishes in the X3 = plane, 



$o(x,i)=/ ^W7V e "'" t+ ' krX| smk 3 x 3 b k + H.c. (2.6) 
/fc 3 >o V 4tHu; 
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where u = |k|. 

Notice that the normalization factor is different for field in half space. 

Writing (ITo(x',i) = $o( x ',^)) and demanding that the equal-time commutation relations are satisfied, 
that is, 

[$ (x,t),n (x / J t)]=^ 3 )(x-x / )-i5 (2) (x|| -xj|)^ 3 + 4). (2-7) 

we can infer that the mode expansion coefficients can be identified with the creation and annihilation oper- 
ators 

[k,bt] = 5 Kk >, [6k,Sk] = 0, [SL,S£] = 0. (2.8) 

Note that the commutation relation above (|2.7p is not in the form normally encountered. For the case of 
a field satisfying the Dirichlet boundary conditions on the x 3 = plane it is necessary to add the "image 
term" to the right hand side such that for points restricted to the X3 < half space where the image term 
has no support the relation (|2.7[) reduces to the standard form. 

With the expression for the free field operator in hand we now calculate the retarded and Hadamard 
propagators. The retarded propagator quantifies the field sourced by the oscillator and describes how energy 
is dissipated from the oscillator into the field 

G% t (x, t; y, H) = iO(t - f )ft>(x, t),& (y, 0] 

= 9(t-t')-^ I d 3 fc-sin(a;(i-t'))e 4k ii- (x ll" y il ) sin(fc3a;3)sin(/s3y3) (2.9) 
2n Jk 3 >o ^ 

/'CO 

= 9(t-t') / dusm(u(t - t')) ■ 7(w;x,y). 



The Hadamard function, which will be encountered in the next section, describes the quantum fluctuations 
of the field and is given formally by the anti-commutator of the field operator: 



G*(x,i; y ,t') = ({*o(x,i),$o(y,i')}> 
1 

~ 2^ 



d 3 k- cos(uj{t - i'))e* k » - (x ii > sin(fc 3 a;3) sm(k 3 y 3 ) (2.10) 
fe 3 >o w 



duj cos(w(i — t')) ■ I(uj; x, y), 



where 

mt/2 f27T 



/(w;x,y) = — 3/ dOj dcj) e lk «< x n-y\0 sin(wcosfe 3 ) si^wcos^). (2.11) 
2"" Jo Jo 



III. EARLY TIME DYNAMICS OF OSCILLATOR-FIELD ENTANGLEMENT 

In Q the dynamics of the entanglement between two inertial oscillators with identical couplings to a 
common massless scalar field was investigated and it exhibited several interesting features. There are roughly 
three temporal regimes if the two oscillators are properly separated and weakly coupled to the field. The first 
regime is at early time, up to t ~ 0(1/ Xq) when accumulated effects of field-mediated mutual influences, 
which manifest themselves in the field's retarded propogators, between the oscillators are still weak and the 
reduced dynamics of the oscillators are dominated by influences of vacuum fluctuations of the field. As the 
effects of field-mediated mutual influences between the oscillators gradually gain strength, the system will 
enter an intermediate regime in which the oscillators' reduced dynamics become quite complicated. In the 
late-time limit, the oscillator approaches a time-stationary state. 

In particular, in Q it was shown that at early time, after causal contact has been established between 
two detectors, an oscillatory pattern of entanglement emerges which varies with the distance between the 
oscillators. This is not due to mutual influences since in the weak coupling limit the mutual influence effect 
is always weak compared to effect of vacuum fluctuation. 
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For our model, we expect a similar division into three temporal regimes. Here in this section we study the 
early time behavior of our model, assuming that initially both the oscillator and the field are in their ground 
states and are uncorrelated. We show that the oscillator-field entanglement also develops an oscillatory 
pattern which varies with the distance between the oscillator and the mirror. At a given instant of time, 
the oscillator-field entanglement exhibits spatial oscillations characterized by frequency Q (we adopt the 
convention c= 1). Additionally, we find that the growth rate at which the oscillator-field become entangled 
oscillates as a function of the oscillator-mirror spacing. 



A. Measure of Oscillator-Field Entanglement 



In our setup, the oscillator and the field together form a closed quantum system which undergoes unitary 
evolution. Since the initial state of the total system is Gaussian and the action is quadratic, the reduced 
quantum state of the oscillator at any time will remain Gaussian. For such systems the amount of bipartite 
entanglement between the oscillator and the field can be quantified by the linear entropy of the oscillator's 
reduced density matrix, which takes on values different from when the oscillator and field are entangled. 
Here the linear entropy is defined as 



S L = 1-V, (3.1) 

where V = Trp\ is the purity of the oscillator's reduced quantum state. Here for our general mixed Gaussian 
state, the purity is given by 

1 1 , , 

V = —== = — (3.2) 

2VdetV 2 J<Q,Q><P,P> - <Q,P> 2 



where < O l5 >^ Tr{p a ■ {O,, 0,})/2, O - (Q, P). 

Mathematically we know that V is always less than or equal to 1. Smaller purity means the reduced 
density matrix of the oscillator is less pure, and correspondingly the oscillator is more entangled with the 
field. 

Another measure of bipartite entanglement between the oscillator and the field is von Neumann entropy, 
which for a single-mode mixed Gaussian state is 

^ = ^ln( — )-ln( — ) (3.3) 

and is a monotonically increasing function of the linear entropy. Both Sl and S v yield the same characteriza- 
tion of mixedness and are equivalent as entanglement measures in the case of a single-mode mixed Gaussian 
state. 



B. Mode Decomposition and EOM for Mode Functions 



We perform the following mode decompositions for Heisenberg operators of the oscillator Q(t) and the 
field with respect to the creation and annihilation operators of the oscillator a, a* and those of the field 
6, 6* : 



<&(*) = A /— [/.(x)a + / a *(x)at] + 



d 3 k 

;3 >o \/2?r3 V 2uj 
d 3 k Ph 



k 3 > 



o V2TT 3 V 2w 



Q+ (i,k)S k + g _(t,k)S^ 

U(x;k)b k + f_(x;k)bt 



(3.4) 
(3.5) 



By plugging (3.4) and (3.5) into (2.2) and (2.3), respectively, we find the following equations of motion 
for the mode functions 
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(8 2 + 2 lQ d t + n?) q a (t) = -^0(t - L)q a {t - L), (3.6) 
(d? + 2 lQ d t + tf r ) q+ (t,k) = -^6(t-L)q + (t-L,k) + ^f 0+ (t,x Q -,k), (3.7) 

in which / 0+ (t,x Q ;k) = e - ia,t+ik ir x ll sinfcgcca and 1q = X 2 Q /(8nM Q ). 

Since we have assumed the oscillator's position to be x<-,(f) = Xg = (0,0, L/2), we have /o+(xg,£;k) = 
e" 8 "* sin^. 

According to our initial condition, the solutions have to satisfy initial conditions that /+(0,x;k) = 
e ik ir x ll smk 3 x 3 , a t /+(0,x;k) = -twe^u*!! smk 3 x 3 , g a (0) = 1, 9 t g a (0) = -ifir, and / a (0,x) = 9 t / a (0,x) = 
g+(0;k)=c> tg+ (0;k) = 0. 

The solution to (|3.6p can be written as the following expansion which is truncated at nL > i, 

9 9 (*) = E4 n) (*), (3-8) 



n=0 



where 



q { °Ht) = qi h Ht), (3.9) 
g£ n >(t) = Jdr 1 G r (t,T 1 )(- 2 ^-)e(t-L)x (3.10) 

drnGrCrn-i - r„)(-^)0(r„ - n£)g< fc >(T n - nL). 



Here gi' 1 '' (i) is the homogeneous solution satisfying the initial conditions 



n=0 



where 



(3.11) 



The solution to (|3.7p can be written in a similar fashion as above, also truncated at nL > t, 

q + (t,k) = Y,q { +\t,-k.), (3.12) 



<4_ 0) (i,k)= f dnGrfrn^gfo+in^k) (3.13) 

^sinM) [( Ml _ M 2 )e— t + (Af 2 e lf2t - Mie" 10 *^*], 



MqVL V 2 

q ( ?- 1) (t,k) = J dr 1 G r (t,T 1 )(-'^)e(T 1 -L)x (3.14) 

dr 2 Gr(ri,T 2 )(-^)fl(T 2 -2L)x 

... / dT n G r (T n -i - T n )6(T n - (n - 1)L)X Q 
Mq/o+(t„ - (n - l)L,x Q ;k), 

where Mi = l/[2(— 00 — ijq + CI)], M 2 = l/[2(— w — i7g — fl)}, G r (t,r) is the retarded Green's function 
which satisfies (<9 t 2 + 2^ Q d t + fl 2 )G r {t, r) = 5(t - r). 
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C. Zeroth-Order Correlation Functions 



As explained before, in order to understand the dynamics of entanglement, it is sufficient to compute 
the covariance matrix of the oscillator's reduced quantum state. For a separable initial state, the covari- 
ance matrix can be decomposed into two parts corresponding to the two sets of operators in the mode 
decomposition: 



< Q(t),Q(t) > = < Q(t),Q(t) > a + < Q(t),Q(t) >, 



fc 3 >o 



(3.15) 



< P(t),P(t) >=< P(t),P(t) > a + < P(t),P(t) > y 



M 2 Q 



fe 3 >0 



(3.16) 



< Q(t),P(t) >=< Q(t),P(t) > a + < Q(t),P(t) > v (3.17) 

= ^liiirnqs) +?:(*)*«.(*)] + J k >o ^^\[ q u^mq + (t,k)+ q+ (t,k)d tq * + (tM- 

Here < Oi,Oj >=< (OiOj + OjOi) > /2. In the weak coupling limit, the effect of reflected influences 
correspond to terms in 13.81 and 13.121 which are of higher than zeroth order, and therefore at early time (up 
to t w 1/7q) accumulated effect of reflected influences is always small. Thus for the purpose of studying 
the early time behavior of entanglement, we can ignore the contribution of reflected influences and restrict 
ourselves to lowest order correlations. 
We have 

<Q(t),Q(t) > a = ak|«£ h) (t)| a 
=^-||(l + fkj£a) e -7gt-«M + i (1 _ n : 4p l)e - 7<J *-Hnt|a j 



<P(t),P(t) >a =^L\d t qi h >(t)\ 2 
|§|i(l + ^)( 7Q + M)e-^-™ + 1(1 - ^)(7 Q - Xl)e-r*««*\\ 



< Q(t),P(t) > a = • d t qi h \t) +qi h) (t) ■ d t q^* (t)) 

= ftfle{[i(l + ^)(7q + <n)e-^*"«« + |(1 - ^)(7q " ffi)e"^+^] 
. [i(i + fk^2a ) e -T«*+«M + i (1 _ fV^raj e - 7<J t-<nt]j_ 

Notice that the zeroth order correlators < ... >a s do not depend on the distance between the oscillator 
and the mirror, but the part induced by vacuum fluctuations does: 



<Q{tlQ(i) > v= 4 >o |i|J-| g W ( t, k )| 2 

= (l&) 2 4^/ rfw l (1 - ^tnr) K A/ i - M 2 )e-^ + {M 2 e int - Af ie - in *) e -^*| 2 . 
Similarly we have 



(3.18) 



8 



< P(t) , P(t) > v = Ml J k3>0 ^\ d t qf (t, k) | 



(3.19) 



(^) 2 iM^!(i- 



| - iu){M x - M 2 )e- Wt + {{in - 7g )M 



2(- 



1q)M 



if 



if2t-j e -7 Q t|2^ 



< Q(t), P{t) > v = M Q f gi ±(qW*( t , k)d t q { °\t 7 k) + q^{t, k)d t gf *(t, k)) (3.20) 



)k 3 >0 2 

\2 1 



fle{[-iw(Mi - M 2 )e- iwt + {{in - 7 Q )M 2 e jm - {-in - 7 Q )Mie- i0t )e- 7 « t ] 
• [(Mi - M 2 )e iult + {M 2 e mt - M^-^e^Q 1 ]}. 

In Q, for the case of two inerial oscillators distance L apart, which are located in free space with identical 
couplings to the field, the zeroth order correlators due to vacuum fluctuation are given as (Mq = 1 and d in 
0| is replaced by L): 



< Q A (t), Q B {t) > v = (4f ) 2 i^ Sduu ^ • 
| (Mi - M 2 )e~ lult + {M 2 e int - Ahe'^e^Q^ 2 , 



< Q A (t),P B (t) > v = (4?) 2 i^ SJ ^TT ■ 

Re{[—iu)(M 1 - M 2 )e- lut + {{in - lQ )M 2 e lUt - {-in - ^^M^-^e'^Q 1 ] • 
[(Mi - M 2 )e iLUt + {M 2 e int - Ahe'^e^Q 1 ]}, 



< P A {t),P B {t) > v = (4f) 2 ^ fduu ^nr ■ 

| - iw(Mi - M 2 )e- lut + {{in - lq )M 2 e int - {-in - 7 Q )Afie- iat )e- 7 «'| 2 . 

Physically, the v-part of the zeroth order correlators < ... > v effectively measure the response of the 
oscillator to vacuum fluctuations of the field. The similarity between the integrands of < Q{t),Q{t) > v in 
Eq. (|3.18j) and < Q A {t),Q B {t) > v in Eq. (|3.21[) is not surprising. In Appendix iBl we show that, if we have 
two inertial oscillators C and D in free space at a distance L apart, with the same coupling constants but 
in opposite signs, then {Qc{t) + Q,o(t))/2 obeys the same equation of motion as the one for Q{t) in our 
model. Thus we see that the self correlator < Q{t), Q{t) > v here has the same value as (< Qc, Qc > + < 
Qd, Qd >)/4+ < Qc, Qd > /2 at t and contains the part of correlations of vacuum fluctuation in free space 
which is odd with respect to the = plane, whereas < Qc{t), Qo{t) > v <x (— Aq)Aq here has exactly the 
same value of — < Q A {t),Q B {t) >v in Q because the two oscillators in Q are identically coupled to the 
field and so < Qa, Qb > v oc XqXq. 



D. Early-time Dynamics of Oscillator-field Entanglement 

With the previous results we can now investigate the evolution of oscillator-field entanglement as the 
distance between the oscillator and the mirror changes. Whereas it may be possible to obtain an approximate 
analytical expression in the weak coupling limit, we can simply study the dependence of the linear entropy 
on L and t numerically, as shown in Figure [TJ 

Our numerical results reveal the following behaviors: 

• For a given L, the entanglement between the oscillator and the field increases monotonically at early 
times, showing that the oscillator is getting more and more entangled with the field after the interaction 
is turned on. (See Figure [1] (upper row) and (lower right).) 



S L vs . L S L vs . t 

t=5 L=l 

Si Sl 




Figure 1. Entanglement dynamics at early time where only zeroth order correlations contribute. Here 7q = 0.02 
and Q — 5. The upper left plot shows the linear entropy as a function of L and t, with the upperright plot being its 
contour plot. The lowerleft plot shows the dependence of linear entropy on L at a given instant of time whereas the 
lower right plot exhibits how the linear entropy evolves with time for oscillator located at a certain distance. 



• Within the light cone, at every given time in this stage, the oscillator-field entanglement exhibits an 
oscillatory behavior with spatial frequency being tt. This shows that the growth rate of Sl with which 
the oscillator gets entangled with the field also oscillates as a function of L, as shown in Figure [T] (upper 
row) and (lower left). 



Such oscillatory behavior in L is solely due to the oscillatory behavior of the zeroth order correlators corre- 
sponding to < Rc, Rd >v; R = Q, P, given in the last paragraph of the previous subsection. After cos 9 in 
&3 = locos 9 has been integrated over the interval [0, 1] in the integration corresponding to < Qc,Qd >v in 
(I3.18P or (13.211) , the field modes on resonance with the oscillator have constructive (destructive) interference 
at the local maxima (minima) of Sl against L at fixed t at early times. 

In Q the early time dynamics of entanglement between the two oscillators exhibit similar oscillatory 
dependence on the distance between them, which is also due to distance-dependent correlations of vacuum 
fluctutaions experienced by the oscillators. More precisely, from (|3 . 2 1 1) one can see that the field modes with 
to = nn I 'L, n = 1, 2, 3, • • ■ , has no contribution at all to the integration of < Qa, Qb >v, and those satisfying 
tanwL = loL give the maximum and minimum values of the factor sinuj L / (u L) in the integrand. In the 
weak coupling limit the integration in (|3.21[) is mainly contributed by the poles at lo w ±0, namely, those 
mode on resonance with the oscillators. So < Qc,Qd > v ~ for L w nir/tt, while < QcQd >v has local 
maximum or minimum values when L is about the solution of tanOi = flL in the weak coupling limit. 
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IV. LATE-TIME STATIONARY LIMIT OF OSCILLATOR-FIELD ENTANGLEMENT 

A. Quantum Langevin Equation and Covariance Matrix at Late-times 

One can compare the dynamics of the oscillator under the influence of the field to the well-studied quantum 
Brownian motion (QBM) model, namely, the retarded and Hadamard Green functions correspond to the 
dissipation and noise kernels respectively and the function I{ui\ x, y) represents the spectral density [see, e.g., 
(22I, [23| ) and the references therein.] In the same vein Eq. (|2.5[) can equivalently be written as a quantum 
Langevin equation 

M Q Q(t) + M Q n 2 Q(t) - jf drn(t, t) ■ Q(t) = i(t), (4.1) 

in which 

,A 

fx(t, s) = X Q 6(t - s) dujI(uj- 1 x Q ,x Q )smuj{t - s), (4.2) 
Jo 

and £(t) = Aq<I>o(xq, t). In the integrals over the frequency above and below we have assumed a high 
frequency finite cutoff A for the quantum field A which regularizes the quantum field's retarded Green's 
function from which we obtain the effective equations of motion of the oscillator In QBM language \i is 
the dissipation kernel, v(t, s) is the noise kernel which quantifies the two-time correlation of the Langevin 
forcing term 

\2 



u(t,a) = ({£(*), £(s)}) =y Q / dw/(w;x Q ,x Q )cos W (f- S ). (4.3) 

Jo 

[From now on we simply denote I(u>) = Xq, Xq).[ Here the average is taken with respect to the initial 
state density matrix of the field. 

By introducing the damping kernel j(t, s) defined by 

d d 

h(t,s) = M q —j(t,s) = Mq—j(t,s), (4.4) 
we can bring f|2 . 5[) to the final form 

M Q $(t) + M Q n 2 r Q(t) + M Q I *ry(jk,r) ■ Q(t) + 2M Q7 (t,t i ) • Q(U) = f(f), (4.5) 



where il 2 , = il 2 — 2MQj(t, t) is the renormalized frequency. 
The general solution to (|4.1[) is given by 



Q(t) = Q (t) + J ds G(t, s)i(s), (4.6) 



where G(t, s) is the retarded Green's function for (|4.ip satisfying 



M Q G{t,s)+M Q WG{t,s)- J dTfj,(t,r)-G(r,s)=8(t-s). (4.7) 

For here and below all the Fourier components are defined for positive frequency only, which correspond 
to Fourier transformation of function defined on t > 0. Accordingly, the Fourier space representation of the 
above Green's function can be written as: 

G(lo) =(M q (-lo 2 + O 2 ) - /iH)- 1 

^(V-^H + n^ 1 !^ 1 . (4.8) 
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In the late-time stationary limit, because /z(t, s) leads to dissipation of the oscillator's free motion, we see 
that 

QH -> G(u)i(u). (4.9) 

As has been explained before, if we assume the initial state of our combined system to be Gaussian, then 
because the total Hamiltonian is quadratic, the quantum state of the oscillator will always be Gaussian and 
can be fully characterized by the covariance matrix 

(vMt) = i<<5i(f)d(t)j + oddity), (4.10) 

in which 6 = (Q,P). 

In the late-time stationary limit the elements Vqp vanish, as can be inferred from Vqp{€) = MVqq{€)/2 
if Vqq(^) approaches an asymptotic constant value. 

The remaining non-zero elements of the covariance matrix at late times are [Tlj : 

V$j = <Q(t),Q(t)> =f A d w G t ( w )./( w ).GH, (4.11) 



t— >oo 

2 l\/r2 r<*t 



Vp x P = < P(t),P(t) > = duj lj 2 M 2 G*(lo) ■ I(u) ■ G(w). (4.12) 

t^oo J 

By applying the fluctuation-dissipation theorem (l8l. [22|. one can eliminate the noise kernel in the covari- 
ance matrix elements, giving 

V QQ = Z [°° <^MG(to)\, (4.13) 



7T 



(I 



Vpp = - duj u 2 MMm[G(uj)}. (4.14) 
^ Jo 

For detailed derivation leading to the above results, please refer to Appendix [A] 

B. Entanglement between Oscillator and Field in Free Space 

We will begin by calculating the late-time oscillator-field entanglement for an oscillator in free space. We 
then compute the oscillator-field entanglement with a mirror present and compare their differences. The 
first step is to regularize the retarded Green's function of the field at the trajectory of the detector. For an 
oscillator interacting with a massless scalar field in free space, the entire system is governed by the action 

S[Q,Q;$,d ll ®] = ^M Q J dt (Q 2 -n 2 Q 2 ) + ^ J dt J d 3 x + X Q J dt Q(t)<S>{x Q ,t). (4.15) 

which is the same as Eq. (|2.ip but without the X3 > restriction. 
Following Q we obtain 

(d 2 + 2 lQ d t + n 2 r ) Q(t) = A Q $(<, x q )/M q , (4.16) 

where Q r is the renormalized natural frequency. 

For free space the retarded Green's function for (|4.1[) is given by 

GH = ^h(^ 2 + Hq? + n?]- 1 , (4-17) 
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where fi r = fi 2 , — 7 2 from which the late-time covariances can be computed: 



V, 



1 



QQ,froo 



dwlm 



^ln 



1 



1 



M g -( w + i 7Q )2 + 02 
7 Q - iVt r 



(4.18) 



PPfrcc 



cLjIm 



Mr 



-(^ + *7q) 2 + ^. 



2?rfi 



Perturbatively in Jq one has 



7q + if2 r 7T 



2A -ml 



7| 

ft 2 



(4.19) 



V, 



QQSrcc 



(4.20) 



PP, free 



Mr 



n r i 



:7q 



q 2 

2(lnA-lnfi r )--£ -1 



which recovers the results in 



(4.21) 



C. Entanglement between Oscillator and Field under Mirror Reflection 

In the presence of a perfect mirror, the entire system is governed by action (|2.1[) . The Heisenberg equation 
of motion for the oscillator after the same regularization as above is 



{d 2 t +2 lQ d t + n 2 r )Q{t)- 



4irL 



0(fi-L)Q(t-L)+\ Q $(t,x Q ), 



and correspondingly 



GH = + *7q) 2 + 0? + (27 Q ^V£)] 



(4.22) 



(4.23) 



where the last term inside the square brackets shows the difference from the free space results. It can be 
interpreted as the contribution from the oscillator's mirror image located at a vertical distance L/2 behind 
the mirror. 

For this case the exact late-time covariance matrix becomes 



V, 



1 



QQ, half — space 



yoo 
PP, half —space 



ttM, 



Q Jo 



dujlm 



1 



Mr 



cfcjlm 



-{u + i lQ Y+^ + {2 lQ e^/L)_ 

uP ' 

-(uj + i lQ )i+& r + {2 lQ e^/L) 



(4.24) 



(4.25) 



Assuming that the oscillator is only weakly coupled to the field, we can perturbatively expand the above 
integrals and get 



QQ, half —space ^QQ, free 



PP,half — space 



PP, free 



(4.26) 
(4.27) 
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where the terms SVqq and 5V£? P represent the first corrections to the covariance matrix elements due to 
the presence of the mirror. Physically keeping only these terms for a single reflection in this perturbative 
expansion is equivalent to ignoring the multiple reflections between the oscillator and the mirror. The exact 
form for the leading order correction is given below: 



ryoo A 



Q Jo 

1 1 



efcjlm 



1 



-2 lo e^ L /L 



7T MnQ r L 



Re 





L \ 







rlj T[0,m r L] 



(4.28) 



5V? P = 



Mr 



dull 



( -2l Q e zujL /L 



M QlQ 

TrQ r L 



Re [(-« + Lfl r )) e mrL T[0, iCl r L]] 



(4.29) 



in the limit of large cutoff A. 

The change of linear entropy due to the presence of the mirror, compared to the case of free space, is given 

as 



AS L = Sl, half —space - S LJrcc == -^Re [e lO " L r[0, iQ r L}] . (4.30) 

In Figure [2] (upper-left) we see that ASl < 0, thus the presence of the mirror always acts to reduce the 
linear entropy between the oscillator and the field, thereby causing them to be less entangled. 

In the same plot one can also see that ASl increases monotonically with L and goes to as L — >• oo. 

This fact can be intuitively understood by inspecting the Heisenberg equation of motion of the oscillator's 
internal degree of freedom. As shown in (|4.22[) , due to mirror reflection, the oscillator's internal degree 
of freedom will have an negative influence upon itself after time L through field propagation. This will 
effectively reduce the quadrature Vqq and Vp*p and hence cause the oscillator to be less entangled with the 
field as it moves closer to the mirror. 

As one would naturally expect, as the oscillator becomes more and more strongly coupled with the field, 
the oscillator-field entanglement increases monotonically, as shown in the lower plot of Figure [5J 



D. Effect of Mirror Image versus Effect of Real Object 



In our setup, the lowest order correction to self-correlators corresponds to the physical process in which the 
oscillator emits a quanta and then interact with it after it is reflected by the mirror. Ths process contributes 
a correction of order O^q/L), according to (|4.28p and (|4.29p . On the other hand, as stated in Appendix 
[Bl in the case of two inertial oscillators with opposite coupling constant, the lowest order correction to self- 
correlators corresponds to the following physical process: oscillator A emits a quanta, which interacts with 
oscillator B, then the back reaction from oscillator B to the field echoes back and interacts with oscillator A. 
The contribution of this process is of order 0(jq/ L 2 ). Therefore, the influence of the mirror and a real image 
oscillator correspond to different physical processes and are of different orders in the coupling constant. 



V. DISCUSSION 



We began with a rather naive question about whether quantum entanglement can exist between a physical 
object (the evil queen Q modeled by a harmonic oscillator) with her mirror image. The immediate answer 
from a formalist could be no, because an image is not a physical object and only physical quantum objects 
can get entangled. What we really wanted to find out in this inquiry is the effect of boundaries on quantum 
entanglement, as that between an atom and its trap or cavity surface. Our calculation leading to an answer 
to this query took on three steps 1) Entanglement exists between an oscillator and a quantum field; 2) The 
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Figure 2. (Upper-left) Leading-order corrections to linear entropy given in (|4.30[) as a function of distance between 
the oscillator and its image due to the presence of the mirror. (Upper-right) Linear entropy of the oscillator as a 
function of distance between the oscillator and its image according to Eqs (|4.24[) and (|4.25[) . (Lower) Linear entropy 
as a function of 7q. Here Mq = 1, 7q = 0.02, Q. r = 5. 

presence of a mirror alters the field configuration, namely, a perfectly reflecting mirror imposes a Dirichlet 
boundary conditions on the field along the mirror surface; 3) The symmetry in the expression we derived 
for the entanglement between the oscillator and the field in half space with Dirichlet boundary conditions 
imposed on the mirror surface suggests that it is as if the oscillator was entangled with its image located 
at distance L/2 on the other side of the mirror. We expound the meaning of these statements with the 
following observations: a) This entanglement is different both from that of an oscillator with the quantum 
field in free space, i.e., without a mirror, and that between two physical oscillators at distance L apart, b) 
The statement of "entanglement between an oscillator with its image" should always be understood as "the 
entanglement of an oscillator with the quantum field in half space under Dirichlet boundary conditions", c) 
Note that this substitute description works only for Dirichlet conditions, it fails if a Neumann condition was 
imposed on the boundary instead. 

Developing this theme further, one can see more clearly which parties are being entangled if one considers 
a microscopic model of a mirror such as the one considered by Galley et al (28| . where the mirror's internal 
degrees of freedom is modeled by an oscillator (called mirosc) with very light mass (whereby the quantum 
fluctuations of the mirrors internal degrees of freedom will be suppressed) . With this setup we could then 
consider three physical degrees of freedom: the oscillator, the field and the mirosc (the light degrees of 
freedom making up the mirror). If one first considers the interaction between the mirosc and the field, this 
would yield in the zero mirosc mass or infinite reflectivity limit the modified field configuration derived here. 
Then from the covariance matrix of the oscillator one can derive the entanglement between the oscillator 
and the modified field. With a microphysical model of the mirror one can calculate how the dynamics of 
the mirosc is altered while interacting with the field, how the field is modified, and how the oscillator is 
entangled with the modified field. The advantage of this is that one can see clearly how this entanglement 
can be interpreted as the entanglement between the oscillator and the mirosc (for an example of the successive 
levels of coarse-graining , see, e.g., (l6|). This reminds us of a similar procedure in electrostatics, namely, how 
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the force between a charge and the induced surface charge density on a conducting plate can be calculated 
using the image charge method. 

Another observation using elementary physics of wave reflection upon a mirror is the following: Since 
the field configuration is at the base of inquires into boundary effects on the oscillator field entanglement, 
the behavior of reflected waves could provide some useful guide in building up our intuition on quantum 
entanglement in this setting. Instead of a mirror with perfect reflectivity we can think of two adjoining 
dielectric media 1, 2 with dielectric coefficients e\ < e 2 (the mirror situation considered above corresponds to 
the case where ei is the vacuum eo, much smaller than e 2 ). For waves propagating from a soft medium 1 to a 
hard medium 2, the reflected wave is inverted. This show up in the reflected field configuration carrying an 
opposite sign from the original field configuration, thus partially canceling it. This cancelation effect is more 
severe near the mirror surface and hence we see the decrease of entanglement as the oscillator gets closer to 
the mirror. If this reasoning is correct then in the reverse situation, if the queen were a water nymph living 
in a lake (medium 2), looking up at the sky (the air is medium 1, where we have assumed e 2 > ei). As waves 
from the heavy medium entering a light medium will be reflected at the interface with a positive amplitude, 
the entanglement would increase as the nymph comes up close to the water-air surface. Thus we can add 
a fourth factor, 4) that of parity in reflection. Fairy tales aside, when the experimental techniques improve 
to the extent that one can measure the quantum entanglement between an atom and the trap surface these 
results could be of some practical value. 

Our next paper (29| will treat the quantum entanglement between an atom and a dielectric medium. We 
will adopt the influence functional method recently used for the treatment of fluctuation forces between an 
atom and a dielectric medium [3. As a small corollary we will be able to check on the correctness of the 
above qualitative argument based on symmetry and parity considerations. Later papers in this series will 
address entanglement domain and entanglement pattern, and a parallel series on quantum entanglement 
in topologically non-trivial spaces starting with R 1 x S 1 [3fJ, which can be applied to atoms in a toroidal trap. 
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Appendix A: Derivation of late-time covariance matrix 

According to our definition, the full-time, exact expression for he QQ-part of late time covariance matrix 



fOO ft ft 

V QQ (t)= / <1uj-I{uj) / d n / dT 2 G( n ) cosKn - t 2 )]G(t 2 ) (Al) 
Jo Jo Jo 

fOO ft /*T2+t 

= / doj-I(oj) / dr 2 / d?G(f - t 2 ) cos[w(f - 2r 2 )]G(r 2 ). (A2) 

JO Jo Jt 2 



At late times the integral can be approximated by 



fOD ft ft 

V QQ (t)~ / du-I(u) / dr 2 / df G(f - r 2 ) cosKf - 2r 2 )]G(r 2 ) 
Jo Jo Jt 2 



(A3) 



pOO pt pT 

/ dio-I{u) df dr 2 G(f -t 2 )cos[u)(t -2t 2 )]G(t 2 ) (A4) 
Jo Jo Jo 

00 ft 

duj-I(u) dTRe{[e- ZUJT G{T)} * [e ltJT G(r)]}. (A5) 



'0 Jo 

When t —> +00, the above approximated expressions become exact, and we have 



V%k = / dw <5» ■ /(«) • G(u), (A6) 
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and similarly, the PP-part of the exact late-time covariance matrix is 



Vp x p = J lo 2 M^G*(lu) ■ I{uj) ■ G(lo). (A7) 
According to previous definitions (|4.2[) (|4.4p we have 



J(w) = -wMqJteWu)], (A8) 

7T 



using the above formula we have 



G*(ui) ■ Re[j(cj)} ■ G{lu) = -(G*(w) • 7*(w) ■ G T (u;) + G*{uj) ■ j{uj) ■ G t (uj)) (A9) 
= 5(^(1 - (-^ 2 + fi*)G»A/ Q ) ■ G T H 
+ Z^-G*(^) ■ (1 - {-u 2 + n 2 r )M Q G T (uj)) 
= ±Im[G(u>)}. 

Therefore we find that, in the late-time limit, we can eliminate explicit reference to the noise kernel and 
express the covariance matrix elements as 

V? Q = f duhm[G(oj)}, (A10) 



V^ P = I dwui 2 -M^lm[G(ui)]. (All) 

Appendix B: Comparison with the case of two inertial oscillators 

In Section VI of Ref. [j| , we have obtained the late-time correlators in the case with two identical Unruh- 
DeWitt detectors at rest at £3 = ±L/2. Let Ao — > — \q for the left detector (Qa) and Ao — ^ +Aq for 
the right detector (Qb), one may wonder whether the detector on the right (Qp at X3 = +L/2) in this 
two-detector case would behave the same as the detector at the same position (say, Qp at x 3 = +L/2) in 
the above single-detector case with its image detector. 

The answer is no. The presence of the other detector separated in a distance L from one detector introduces 
corrections to the late-time correlators of a single detector, which are O^q/L 2 ) for the self correlators and 
O^q/L) for the cross correlators. This is different from the above self correlators Vqq and Vpp, which have 
SVqq and 5V PP in 0(j Q /L). 

This can be understood as follows. The late-time behavior of the correlators in Q are determined by mode 
functions (i, k) and q B (t, k), whose equation of motion reads (Eq.(13) in Q with d and Aq modified) 



(O 2 + 2 lQ d t + n 2 r ) q ( B + \t,k) = —^9(t - L)q£\t -L,k) + \ Qe -i"t+ik 3 L/ 2j (B1) 
{d 2 t + 2 lQ d t + n 2 r ) q { +\t, k) = -^6(t - L)q { +\t -L,k)- X Q e-^- ik ^ 2 . (B2) 

Similar to Eqs. (A6)-(A11) in Q, these equations give the same c£ for q^ like (^4.11), while the counterpart 
for Qa is — c£. So the late-time self correlators are still given by Eqs. (48) and (50) in @, and the cross 
correlators are those in Eqs. (49) and (51) multiplied by — 1. From Eqs. (54) and (55) in Q, one can see that 
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in weak coupling limit the late-time cross correlators are 0( 7 q/L) and the correction to the self correlators 

is 0( 7 |/i 2 )- 

On the other hand, the detector in this paper has 



(df + 2 lQ d t + tf r ) q + (t, k) = -^0(t - L)q + (t - L, k) + ^-e"-' sin 



where q+{t — L,k) on the right hand side can be interpreted as the image of q+. Indeed, (J2HJ gives 



(B3) 



9+(t,k)| 



sm 



k,L 



/M Q 



t»i/ 7Q w 2 + 2i lqUJ _ ^2 _ (2 7g e iwi /L) ' 



so at late times, 



< Q(t),Q(t) > v 



d 3 k 



(2tt) 3 2w 
ft 



l<z+(*.k) 



rfw Im 



2* 7qW - n» - (2 7Qe -*-i/L) 



(B4) 



(B5) 



which is exactly the Vqq in (|4.24p after letting h = 1 . 

Note that the first order correction to Vqq in free space has exactly the same value as the late-time cross 

correlator ( {Qa,Qb} ) in the case with two inertial detectors considered above. However, the latter will 
not enter the reduced density matrix of Qb- Compare (|Z?3|) with (|i?l[) and (ji?2[) . one can see that it is 
(q^ + Qjl )/(2i) rather than has the same late-time behavior as q+ for Mq = 1. 
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